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f^ , A theory is developed of the intricately fingered patterns of flux domains 
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observed in the intermediate state of thin type-I superconductors. The pat- 
terns are shown to arise from the competition between the long-range Biot- 
Savart interactions of the Meissner currents encircling each region and the 



'^ ' superconducting-normal surface energy. The energy of a set of such domains 



is expressed as a nonlocal functional of the positions of their boundaries, and a 
simple gradient flow in configuration space yields branched flux domains qual- 



H ' itatively like those seen in experiment. Connections with pattern formation 

in amphiphilic monolayers and magnetic fluids are emphasized. 
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When a thin film of a type-I superconductor is placed in a magnetic field normal to 
the sample, the large demagnetizing effects associated with the film geometry preclude the 
establishment of the Meissner phase (with magnetic induction B = 0). The sample instead 
accommodates the field by breaking up into a large number of superconducting (B = 0) 
and normal (B 7^ 0) regions, usually forming very intricate patterns []I|. These arise in part 
from the competition between the magnetic field energy of the domains and the surface 
energy between the superconducting and normal regions. However, they have not been 
successfully explained theoretically, for virtually all theories of these patterns, starting from 
that of Landau [Q, have explored this competition with variational calculations that assume 
reg'«/ar geometries of the fiux domains [0]. The hypothesized parallel stripes, ordered arrays 
of circles, etc. are rarely seen, the norm instead being the disordered patterns documented 
in classical experiments [0,^,0. Moreover, the temperature and magnetic field history of 
the sample strongly infiuence the patterns, suggesting that they are likely not in a global 
energetic minimum. 

Recent work in this area has emphasized that the diffusion of magnetic flux in the normal 
phase can influence the domain morphology 0,0. From asymptotic methods applied to the 
time-dependent Ginzburg-Landau (TDGL) model [§|,3, a free-boundary representation for 
the motion of superconductor-normal (S-N) interfaces has been derived which is nearly 
identical to those describing the growth of a solid into a supercooled liquid, where the 
interface motion is unstable (producing dendrites, for example). By analogy it was suggested, 
and conflrmed by numerical solution of the TDGL equations f^J^, that the growth of the 
superconducting phase into the supercooled normal phase should be dynamically unstable, 
leading to highly ramifled domain shapes. While such diffusive instabilities may play a role 
in the pattern formation in the intermediate state, these studies are not directly applicable 
as they have all ignored demagnetizing effects. 

Here we take a different approach to the intermediate state, ignoring entirely the diffusion 
of flux and focusing instead on the demagnetizing flelds. We address the basic question: 
What is the energy of a thin multiply- connected superconducting domain, the normal regions 



of which are threaded with a magnetic field? A central issue is the degree to which the 
interactions between the Meissner currents, flowing along the S-N interfaces within the film, 
are screened by the superconducting regions. Pearl [|iy] and Fetter and Hohenberg ||rT| made 



the important observation that vortices in a thin film interact with a potential V{r) ~ 1/r for 
large separations r, while for small r, V{r) ~ ln(A/r), with A an appropriate cut-off. Unlike 
in bulk, the interactions are unscreened, being dominated by the electromagnetic fields above 
and below the film. This suggests the simple model developed here: domains bounded by 
current loops interacting as in free space, endowed with line tension, and subject to the 
constraint of constant total magnetic fiux through the sample. We focus on how such loops 
relax toward a local energy minimum when prepared in a nonequilibrium configuration. Such 
a model allows for the study of the shape of the domains, but does not readily account for 
changes in topology (e.g., fission). We conclude that the long-range interactions destabilize 
fiux domains of regular shape, producing branched, fingered structures as seen in experiment. 
Apart from the global fiux constraint, this model is equivalent to one for domains of 
magnetic fiuids |T2| in Hele-Shaw fiow ]IB| , P^ , where experiments ||T3|,|TB[ show patterns like 



those observed in the intermediate state, with the same history-dependence noted earlier. 
Thin magnetic films |T^ and monolayers of dipolar molecules |[TB| exhibit similar behavior. 



and are described by such models through the underlying similarity between electric and 
magnetic dipolar phenomena [|l9| , p0| . This model is also very similar to one introduced for 



pattern formation in certain reaction-diffusion systems [21]. There, fronts between regions of 



different chemical composition move in response to line tension and a nonlocal Biot-Savart 



coupling, and may exhibit labyrinthine patterns, as observed in recent experiments p2 . 

To begin, we note a crucial separation of length scales between the typical size of the 
flux domains (ca. 0.1 mm) and the penetration depth (ca. < 1 fim) [|T|. Thus on the scale 
of the patterns the superconductor-normal interface is sharp, and we may view the order 
parameter magnitude as piecewise constant. It follows that the energy is determined solely 
by the locations of the S-N interfaces. 

In this macroscopic approach, the energy JF of a conflguration of flux domains arises 



from the bulk free energy of the two coexisting phases, the boundaries between them, and 
the magnetic field energy, 

^ = J^hvLlk + ^int + Afield ■ (1) 

Suppose the film has thickness d, total area A, volume V = Ad, and contains a set of normal 
domains i with area Ai, length Li, and whose boundary positions are rj(s) (Fig. ^. We 
assume r^ is independent of z, neglecting the "fanning out" of the domains near the film 
surfaces @. The two phases occupy volumes Vs and Vn = dY^^Ai, with Vs + Vn = V. Their 
bulk free energy densities Fg and F„ define the critical field Hc{T) as F„ — Fg = H^/8tt. 
With pn = An/ A the area fraction of the normal phase, o"sn = {Hl/9>Ti)^ the S-N interfacial 
tension (where A(T) can be interpreted as the interfacial width), we have 
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Ad- J-.,, = :|^ A ^ L, . (2) 



The complexity of this problem lies entirely in the computation of the field energy. An 
applied field H = HaG^, produces a field Hn^z in the normal regions, where AH a = AnH^ 
by flux conservation, so 

Hn=—. (3) 

Pn 

The requirement of tangential continuity of H across a S-N interface gives the field in the 
superconducting region as H^ = -ffgC^ = Hn^z. Taking the superconducting regions to 
be perfectly diamagnetic (with susceptibility x = — l/47r) yields a magnetization M = 
— (if„/47r)e2, which guarantees that the magnetic induction in the superconducting regions 
vanishes: B^ = H^ + 47rM = 0. The field energy is that of domains of total magnetization 
Mtot (as yet unknown) in an external field H^, 

Afield = -- /rf'VMtot ■ H, . (4) 

We quickly point out that if we ignore the fringing fields entirely and assume Mtot is 
completely uniform, the field energy would be 



Minimizing the sum J-'huik + -Afield s^'^ with respect to p„, we obtain the standard result 
PB| for the equihbrium value p* = Ha/ He] i.e., if„ = ifc- But we obtain no information 
about the shapes of the domains. The inclusion of line tension alters slightly the equilibrium 
area fraction, and renders the equilibrium domain shape a circle, as that has the minimum 
perimeter for a given area. 

A field energy sensitive to the shape of the flux domains must account for fringing 
magnetic fields. Since the calculation of these is well-documented [13,14], we summarize 
only the essential results. First write the total magnetization as Mtot = xHcap, where 
Heap is the magnetic field arising from the uniform magnetization x^n- The subscript 
"cap" indicates the equivalence of this problem to that of a uniformly charged parallel plate 
capacitor. Now write H^ap = H*^*^) + H'^^\ where H'^°^ = H^^^e^ is uniform between the plates 
and vanishes elsewhere. The choice H^'^' = —AirxHn yields V ■ H*-^-* = 0, so H*^^-* = V x A, 
for some vector potential A. Choosing an appropriate gauge, we find V^A = V x H*^°\ or 

A(r)^-§i:/>.f^^^^|^. (6) 

where tj(s) is the tangent to rj(s). Substituting into (H), performing the integrations over 



the thickness of the slab |14|, the energy of a set of interacting current loops is found to be 
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Y, jds jds ti ■ ij ^ij \ , (7) 



where h = Ha/H^ is the reduced applied field, tj = tj(s), and $jj = (1 + 6ij)^{Rij/d), with 
Rij = |rj(s) — rj{s')\. The interaction potential 



$(0 = sinh-l(l/0+e-Vl + e^ (8) 

is Coulombic for ^ ^ 1, $ ^ 1/(2^), reflecting the usual self-induction energy for elementary 
currents, while for ^ ^ 1, $ ?^ ln(2e^^/^), with the film thickness h acting as a cutoff (as in 
the vortex potential V{r) above). 



The motion of S-N interfaces is determined by the normal component of the force acting 
at a point rj(s) on a given interface. It is the negative functional derivative of the energy 
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where Qi = daVi-doYi is the metric, /Cj(s) the curvature, lij the unit normal, Rjj = R,;i7|Ri 
$'(.^) = 1 — (1 + 'C^^)^^^, and n is an effective pressure given by 



n 



p2 
rn 






1 



(10) 



Equation @ reveals that the pattern formation arises from a competition between a 
magnetic pressure which incorporates flux conservation, the Young-Laplace force due to 
interfacial tension and the Biot-Savart force from the circulating currents. This additional 
long-range contribution to the boundary force has been noted before for magnetic fluids 



|T3|,li| and amphiphilic monolayers [Q. From extensive studies [13,n,15,16,20| of this 
competition in domains of fixed area, we know of the existence of branching instabilities of 
circular domains and buckling instabilities of stripes, phenomena which should carry over 
to the present problem with nonconserved area. Indeed, buckled domains are well known in 
type-I superconductors [H,^ and are seen as well in amphiphilic monolayers [^ . 



The simplest dynamical model for interface motion balances the normal component of 
the force against a viscous force localized on the interface, yielding an equation of motion 
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where r^ is a kinetic coefficient. This equation of motion is local in time; in the present context 
it arises from the neglect of the diffusion of the magnetic flux in the normal phase , which 
is equivalent to assuming that the normal state conductivity is zero. The magnetic vector 



potential is then "slaved" to the order parameter, a limit considered also in the reaction- 
diffusion context 1^. To estimate rj, we appeal to earlier results in the slaving limit for 
strongly type-I bulk superconductors and obtain 
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where T^ is the critical temperature and ^o is the zero-temperature correlation length [[26 
These results can be used to estimate the time scales for domain motion from other experi- 
mentally measured parameters. 

The dynamical evolution of flux domain shapes governed by Eqs. (|9|)- ([TI|) is a many- 
body problem of considerable computational complexity. To gain insight into its behavior, 
we consider here a mean-field description in which only the self-interaction of the current 
loops is considered; the Biot-Savart interactions between loops are neglected, while the 
flux-conservation condition is enforced by assigning the loop to a cell of area Aceii, with 
Pn = An/Aceii- By appropriate rescaling of the spatial variables, we deduce from Eqs. (^ 
and ( pUD that the location of the system in the H — T plane is uniquely specified by the two 
dimensionless quantities h and A/AJ^-^. 

To see the effects of the Biot-Savart interaction on the stability of a circular flux domain. 
Figures §a and b show the evolution of a domain prepared with an area significantly less than 
the equilibrium value. In Fig. ^a the Biot-Savart coupling is omitted, and the circle simply 
relaxes to a new radius driven primarily by the magnetic pressure. Including the long-range 
interactions (Fig. 0b) results instead in the formation of a branched flux domain with three- 
fold coordinated nodes, very similar to those seen in experiment [||,H,0. In the early epoch 
the shape evolution is primarily a dilation with little change in shape, while the branching 



instabilities occur on a longer time scale. As found in previous studies ||T6[, in the later 
stages of the shape evolution the driving force for the interfacial motion becomes extremely 
small, with very small energy differences between rather different shapes. This suggests 
that the interface motion would be extremely sensitive to external perturbations such as 
impurities or grain boundaries, which would then be effective in pinning the interfaces, not 



unlike the pinning of vortices in type-II superconductors ||2^. This may contribute to the 
history- dependence of the patterns discussed in the introduction. 

Elsewhere we will discuss the derivation of the present model from the more fundamen- 
tal TDGL equations by generalizations of existing asymptotic methods [|,^^, as well as 
the details of stability analyses for various regular geometries. Given the connections out- 
lined here between pattern formation in the intermediate state and in other systems such as 
amphiphilic monolayers and magnetic fluids, we suggest that it is of interest to extend exper- 
imental studies of flux domain shapes to probe systematically [|l^ the branching instabilities 
as a function of applied field and temperature. 
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FIGURES 
FIG. 1. A thin slab of a type-I superconductor, of thickness d, viewed along the applied field 

Hft. The normal regions are shown shaded. Adapted from [Qj. 

FIG. 2. Numerical results from the gradient-flow model for the evolution of a flux domain 
boundary. Dashed initial condition is a circle of unit radius perturbed by low-order modes, evolving 
with A = 0.01 and h = 0.5. (a) Rapid relaxation to a circle in the absence of long-range forces, (b) 
Rapid dilation followed by fingering on a much longer time scale in the presence of the Biot-Savart 
interactions. 
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